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Abstract 



The infinite-dimensional spin one-half Falicov-Kimball model in an external 
magnetic field is solved exactly. We calculate the magnetic susceptibility in 
zero field, and the magnetization as a function of the field strength. The 
model shows an anomalous magnetic response from thermally excited local 
moments that disappear as the temperature is lowered. We describe possible 
real materials that may exhibit this kind of anomalous behavior. 
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I. INTRODUCTION 



The spin-one-half Falicov-Kimball mo defi was introduced in 1969 to describe metal- 
insulator transitions in transition-metal and rare-earth compounds. The metal-insulator 
transition is a charge-transfer transition, where the system consists of both localized (insu- 
lating) electronic states and delocalized (conducting) electronic states, and the transition 
occurs when the electron filling switches from the localized to delocalized states (or vice 
versa). This is the simplest possible metal- insulator transition, since the character of the 
individual electronic states does not change at the transition, rather it is just the occupancy 
of those states that varies. 

The localized electrons are strongly interacting with each other, so double occupancy 
of the localized orbitals is forbidden. The conduction electrons, however, are chosen to be 
noninteracting, since the Coulomb interaction between them is screened. The only remaining 
interaction term in the Falicov-Kimball model is the mutual repulsion between a conduction 
and localized electron that sit on the same lattice site. For any fixed configuration of the 
localized electrons, the quantum-mechanical problem for the conduction electrons can be 
solved in a single-particle basis. The many-body aspects enter from taking an annealed 
thermal average over all possible configurations of the localized electrons. 

We are most interested in the situation where the localized level lies just above the Fermi 
energy, so there are no localized electrons at T = 0. As the temperature rises, then it becomes 
more favorable to thermally occupy the localized levels because of the gain in entropy of the 
local moments. This will have a dramatic effect on the uniform spin susceptibility, as the 
localized electrons will generate a Curie-like response, which is much larger than the Pauli- 
like response of the conduction electrons. The net effect is that the susceptibility will show 
a peak at a characteristic temperature which is determined by the point where the localized 
electron occupancy saturates. In addition, we expect there to be interesting behavior as 
a function of magnetic field, where the proximity of the localized levels close to the Fermi 
energy can induce a rapid change in the magnetization, producing metamagnetism. 

The Falicov-Kimball model was initially applied to a wide variety of transition-metal 
and rare-earth compounds such as SmB 6 , V 2 3 , NiS, etc. which display both continuous 
and discontinuous metal-insulator transitions as a function of either temperature or pressure. 
However, the Falicov-Kimball model has not been generally accepted as explaining the metal- 
insulator transitions of all of these different materials for two reasons. First, it neglects all 
effects of hybridization, and hence all Kondo-effect physics, and second it does not take into 
account the fact that most of these materials also undergo a structural phase transition. 

Recently, however, a new material was discovered that undergoes a classic charge-transfer 
metal-insulator transition without any structural phase transition — it is Nil 2 when put under 
high pressure^. This layered compound crystallizes in the CdCl2 crystal structure, in which 
the Ni ions form a close-packed plane, with close-packed I planes lying both above and 
below. The Nil 2 planes are then stacked in an alternating packing scheme. The Fermi level 
of Nil 2 lies within the Ni d-bands, since the Ni donates one electron to each I ion that lies in 
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the plane above and below. The ground state is an antiferromagnetic insulator. As pressure 
is applied, the I p-bands move closer to the Fermi level causing the Neel temperature to 
rise, until they reach the Fermi level, and electrons from the I ions move onto the Ni ++ ions 
changing them to Ni + and quenching the local magnetic moment. The remaining holes in 
the I p-bands are conducting, and the material becomes metallic. 

Another class of materials that might be described by an effective Falicov- 
Kimball model is the class of Yb-based valence-fluctuating (VF) compounds (like 
YblnCu4, YbIni_ rit Ag :r Cu4 and Ybi_ a; Y a; InCu4) which exhibit large anomalies in their 
thermodynamic! 1 !, spectroscopic!!!^ and transport properties!'!!, around a characteristic, 
sample-dependent temperature Ty (40 — 60K). The magnetic susceptibility xCO °f these 
compounds is of a Curie- Weiss form at high temperatures, and the anomaly appears!'!!, 
as a pronounced asymmetric peak just above this characteristic sample-dependent temper- 
ature Ty. Below Ty, x{T) exhibits Pauli-like behavior, while at the lowest temperatures 
a weak Curie-like upturn is often seen. The X-ray analysis shows that the Yb-based VF 
compounds crystallize in a C15b structure, that the lattice constant a(T) changes at Ty by 
a small amount!iE, and that there are no structural changes. There is no evidenceS! for 
magnetic ordering below Ty, instead the Yb ions fluctuate between a 2+ and 3+ state, but 
the average volume change at Ty is much less than what one would expect from a com- 
plete Yb 3+ -Yb 2+ transition.! The anomaly in the high-field magnetization of YbInCu4@ and 
YbIn 1 _ x Ag 1 .Cu 4 appears as a sudden increase of the slope and the saturation of M(H) at 
about Hy ~ 30 — 50 T. The Zeeman energy at Hy (at T <C Ty) is comparable to the ther- 
mal energy, at Ty (at H <C Hy). At elevated temperatures, this metamagnetic transition 
disappears!. The magnetostriction data! indicate that the metamagnetism relates to the 
valence fluctuations of the Yb ions but the field-induced change of the /-electron valence is 
even smaller than the temperature-induced one. 

The large magnetic anomalies accompanied by small lattice changes can be reconciled 
by a model in which there is some intrinsic disorder between Yb and In sites, such that 
a minority of Yb ions (in the C15b lattice) lie at "ill-placed sites" . The correctly placed 
Yb ions are hybridized with the ligands and are in the VF state at all temperatures. The 
"ill-placed" Yb ions are unhybridized (because the Yb-In distance of an "ill-placed" Yb ion 
[that sits on an In site] is nearly twice as large as the Yb-In distance for an Yb ion that sits 
on the correct lattice site) and undergo a charge-transfer transition at Ty (one 4f electron 
jumps into conduction band and leaves behind a magnetic hole). As a consequence, the 
susceptibility acquires a large asymmetric peak just above T v , while the low-temperature 
magnetization exhibits a metamagnetic transition for magnetic fields of the order of ksTy. 

This recent experimental work provides our motivation to study the magnetic response 
of the spin-one-half Falicov-Kimball model, and to compare the results with the available 
susceptibility data for NH2 and YbInCu4. On the theoretical side, the Falicov-Kimball model 
has also received renewed interest, with work concentrating on the linear and nonlinear char- 
acteristics of the modelB and on the intermediate- valence and metal-insulator transitions!!!. 

In Section II we present the formalism, showing how to determine the Green's functions 
in a magnetic field, and the different magnetic and charge- density- wave response functions 
in zero field, and the numerical methods employed in the exact solution of the problem. 
Section III presents our results for the infinite-dimensional hypercubic lattice and Section 
IV summarizes the possible applications to real materials. Our conclusions appear in Section 
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V. 



II. FORMALISM 

The Hamiltonian of the Falicov-Kimball model consists of two types of spin-one-half 
electrons: conduction electrons (created or destroyed at site i by d\ a or d ia ) and localized 
electrons (created or destroyed at site i by fj c or f ia ). The conduction electrons can hop 
between nearest-neighbor sites, with a hopping matrix —t^ =: — t*/2y/D, where we have 
chosen to examine hypercubic lattices in D-dimensions, and we choose a scaling of the 
hopping matrix that yields a nontrivial limit in infinite-dimensions (D — > oo )0. The f- 
electrons have a site energy Ef (which will be chosen to lie just above the Fermi energy at 
T = 0), and a chemical potential /i is employed to conserve the total number of electrons 
n d] + n di + n /t + n fi = constant. Finally, there is a restriction of no double occupancy of 
the /-electrons on any site (implying the Coulomb repulsion Uff between two /-electrons is 
infinite) and there is a Coulomb interaction U between the d- and /-electrons that occupy 
the same lattice site. In addition, we include a coupling to an external magnetic field h with 
a Lande g-factor added for the localized electrons. The resulting Hamiltonian i 

i,aa' i,<J 

-VBhJ2°( 2 4*dia + 9 fl fia)- (1) 

We solve the Falicov-Kimball model by working in the infinite-dimensional limit. The 
bare density of states p(e) on a hypercubic lattice becomes 

p(6) = ^exp[-e 2 /r 2 ], (2) 

and we take t* as the unit of energy (t* = 1). (The alternative lattice to consider is the Bethe 
lattice, where the bare density of states becomes Wigner's semicircle. The latter density of 
states is closer to a three-dimensional band because it increases like y/E at the band edges. 
Since the Fermi level always lies far away from the band edge in all calculations considered 
here, we don't expect there to be much difference between the solution on a Bethe lattice 
versus on a hypercubic lattice. In fact, we have checked that this is not the case for our 
results in Section III.) In the infinite-dimensional limit the local approximation becomes 
exact0 so that the self energy has no momentum dependence, and the local conduction- 
electron Green's function satisfies 

G a {iuj n ) =: G na = [ de, (3) 

where u n =: 7rT(2n+l) is the Fermionic Matsubara frequency, and p(e) is the noninteracting 
density of states. The self-energy is determined by a three-step process: (i) first, the site- 
excluded Green's function G®.(iu n ), which holds the information about all of the other lattice 
sites, is found by adding back the self energy, 
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G ° ( ^ = <h + eT ; (4) 



-'no- 



then (ii) the local Green's function is calculated from a weighted sum over the different 
possible localized-electron states, 

T I 

where Wq = 1 — w\ — w\ and wl are the localized electron occupancies for no f-electrons 
(wq) and an f-electron with spin a (wf); and finally (hi) the self-energy is determined by 
subtracting off the "bare" Green's function 

= [G^(iw n )] -1 - G~l. (6) 

The f-electron occupancies are found by explicitly calculating the path-integral for the parti- 
tion functional (generalized to include the magnetic field) and are determined by w = Z /Z, 
wl = Zl/Z, and Z = Z + Z[ + Z{, with 

oo 

z = IJ z\z^ 

n=—oo 

1 oo 

Z{ = exp[-P{E f + U-ii- -gfiBcrh)} J[ {Z\ - U){Zi - U), (7) 

n=— oo 

where Z a n := [Gl{iu n )\-\ 

The algorithm for finding the self-consistent solution to these equations is the same as 
was used in the Hubbard mo delS: (i) begin with the self energy set equal to zero; determine 
the local Green's function from Eq. (|3]); (iii) determine the site-excluded Green's function 
from Eq. (f|) and the f-electron occupancies from Eq. (^); (iv) determine the new local 
Green's function from Eq. (||) and the new self energy from Eq. (^j. This new self energy 
is plugged back into step (ii) and the process is iterated until it converges (which typically 
takes less than 100 iterations for convergence to one part in 10 6 ). 

In addition to finding the Green's functions in a magnetic field, we can also determine 
the susceptibilities to both charge- density- wave order and spin-density-wave order in zero 
magnetic field. The derivation of these quantities is similar to that of Brandt and MielschEl, 
and we simply include details not presented before and summarize the generalizations needed 
when the electrons have spin. 

We begin by adding an ordering field Y2j hjdj a dj a to the Hamiltonian and examining 
the linear response in the limit as h a - — > 0. The real-space conduction-electron correlation 
functions can be expressed in terms of the Green's functions via 

(K - - «,»> = — T 2 £ — 2*. (8) 

n an 'j 

Following the standard techniques^ one can express the correlation functions in terms of 
charge- density- wave (CDW) (%) and spin-density-wave (SDW) (%) susceptibilities in mo- 
mentum space 
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m 

X dd {q) = X dd \q) -TY,X d n d0 (q)r d n d m X%(q), (9) 

m 

where the susceptibility is found by summing over Matsubara frequencies [x dd {q) — 
TJ2 n xffiq)}, xff°(q) 1S the bare particle-hole susceptibility, 



X?f(q) = -Tj2G n (k + q)G n (k), 



k 

oo 



dy 



p(y) 



iuj n + n - S n - y 



dz — 

00 lUJr 



fi-E n -X(q)y-zJl-X^q) 



(10) 



where all of the wavevector dependence is included in the term X(q) = 2j=i( c °s qj) / d (the 
Green's function and self energy in zero magnetic field have no spin dependence), and T 
i^nL) are the irreducible CDW (SDW) vertex functions, 



M 
nm 



-\dd 



1 

T 



dG 



+ 



m\ 



dG 



1 

T 



n t 



dG 



dG. 



mj. _ 



(11) 



Now the self energy can be expressed as an explicit function of the Green's function (of the 
same Matsubara frequency) and the f-electron filling [by substituting Eq. (§) into Eq. @ 
and solving the resulting quadratic equation for £ n ] 



2G r , 



^±4 

2 2,G r 



1 - 2(1 - 2w 1 )UG ntT + £/ 2 G^ 



(12) 



t I 

with W\ = w{ + wj; the total f-electron concentration. The difficulty in calculating the 
derivatives (to find the irreducible vertices) arises from the fact that the f-electron filling w\ 
is an explicit function of all of the Green's functions. 
The CDW vertex becomes 



<9£ 



T \dG n1 



&mn 



1 /<9£ nT 



T V dw! 



G 



dwi dwi 
+ 



dG 



m]~ 



dG 



ml 



(13) 



In the zero-magnetic-field case, the two derivatives in the square brackets are equal to each 
other, which simplifies the analysis. Substituting into the Dyson equation Eq. @ then yields 



Xt d (q) = XT (?) 



ddO, 



1 + xT 

with the function 7(g) defined by 

7(9) =: £ Xf(l) 



(14) 



<9u>i 9wi 

+ 



<9G 



ml 



dG 



ml 



(15) 



Multiplying Eq. fll4|) by dwi/dG n ^ + dwi/dG n i and summing over n yields an equation for 
7(g). Brandt and Mielsch show how to massage that formula into a form where the individual 
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derivatives can be explicitly calculated. The algebra is somewhat long, and entirely contained 
in their work. The final result for 7(g) is 



/ dwi , dwi \ [ ur a ( as n \ 

1 i 7~t „ r'\ n2 I £)£„ \ 



7(5) 



l+G nVn (q)-G 2 J§^ 



,)Z 



l+G nVn (q)-Gl( 



(16) 



with r/(g) defined by 



r) n {q) =: G r , 



1 



X d n d0 (q) 



(17) 



Finally, the CDW susceptibility is computed from 



X M (q) = -T£ 



n 1 + G n T] n (q) - G 2 n 



Ull 



(18) 



This susceptibility diverges when 7(g) diverges, which occurs when the denominator of 
Eq. (jlB]) vanishes. This yields the same result as Brandt and Mielsch for the spinless case, 
except for an additional factor of 2 multiplying the sums in the numerator and the denomi- 
nator arising from the two derivatives of w\ (which are both equal). 

Each of the derivatives appearing in Eqs. ([Ll]) and ([TJ) can be directly calculated. A 
straightforward differentiation and simplification [employing the quadratic equation whose 
solution gave Eq. (|T2|)1 yields 



dwi dw\ 



+ 



2w 1 (l-w 1 )UGl 



;i + G* n E n )(l + G n [£ n -f/])' 



(19) 



C 2 (8Kn\ (l + GA)(l + G»[E„-[/]) 
n \dG n ) wi l + G n (2E n -U) ' 1 > 



and 

2 / (9S n \ UG 2 n 

n \dwJ Gn l + G n (2E n -UY 1 ) 

For the cases examined here, we never found any divergences of the CDW susceptibility. 

We must follow a similar procedure to find the SDW susceptibility, but the algebra 
simplifies tremendously in zero external field, because the irreducible SDW vertex becomes 

T dd _}_( 9Z n A 



T \dG nh 



u-'i 



since the derivatives with respect to Wi are equal and hence cancel. The SDW susceptibility 
then assumes the simple form 
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X M {q) = ~T £ - - : fBE x , (23) 



Wl 



which is never expected to diverge at any finite temperature. 

Brandt and Mielsch also show how to calculate the mixed susceptibilities that correlate 
the d-electron charge (or spin) with the f-electron charge (or spin). The idea is that deriva- 
tives of with respect to the ordering field h a - produce the mixed susceptibilities, so a 
direct calculation of the cM-susceptibility involves total derivatives of the self energy with 
respect to the Green's functions, which can be expressed in terms of the irreducible vertex 
(as done above) to find the susceptibility, or can be expressed as partial derivatives with 
respect to the Green's function and the /-electron filling (since the self energy is a function 
of G n and w\). Hence one can solve directly for the mixed susceptibilities (the details appear 
in Brandt and Mielsch's work). The final result for the mixed CDW susceptibility is 



X df (q) 



1 + G n r] n (q) - G 2 r 


( d£„\ 
\dG n J Wl _ 


X d n d {q) + Gl 


G 2 J 







7(9), (24) 



in which the Matsubara-frequency dependence actually drops out of the expression [after 
using Eq. fllj])], and it is equal to 7(g). Similarly, we find the mixed SDW susceptibility 
vanishes x (q) = 0- 

Now, we could have also calculated the mixed susceptibility by adding an ordering field 
for the /-electrons, and differentiating the (i-electron Green's function with respect to the 
/-electron ordering field. The results of each of these calculations must be the same, since 
the <i/-susceptibility is equal to the /d-susceptibility. However, if we write out the expression 
for the /^-susceptibility explicitly in terms of the derivatives with respect to the /-electron- 
ordering field, we find that it includes a dependence on the //-susceptibility. Hence the 
//-susceptibility can be determined from the (i/-susceptibility found above. In the CDW 
channel, one finds that the //-susceptibility satisfies 



G'n 

E 



1+G n r, n (q)-Gl(^) w 



(25) 



This result can be further simplified following the same steps as above, to rewrite the deriva- 
tives in terms of explicitly calculable quantities. After a significant amount of algebra, similar 
to what was done for the <i<i-susceptibility in the CDW channel, we find 

vffta) = (26) 



Er 



l+G nVn (q)-Gl(^ 



which has the same denominator as the dd- and <i/-susceptibilities. The CDW transition 
temperature, is thereby determined by the condition that this denominator vanishes. 

A similar analysis for the SDW channel does not yield any useful information because 
X d *{q) — 0. Instead, we can directly compute the uniform (q = 0) //-susceptibility for SDW 
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order by differentiating the fillings w° with respect to the /-levels EJ. The calculation is 
simplified by the fact that the factors Z° depend only on the total /-electron filling w\ which 
does not change in a magnetic field (only the relative fillings w\ — Wi change), and the final 
result is 

X"(0) = ^, (27) 

which is in the Curie- Weiss form. 



III. RESULTS 

We present the magnetic response of our numerical solutions for a variety of different 
cases in Figure 1. Each figure corresponds to a different value of total electron concentration 
(n = rid + nf), with ten values of Ef/t* ranging from to 4.5 in steps of 0.5. Fig. 1(a) is 
a typical high-density result. There are 2.5 electrons per impurity site, implying nf ^ 
at all temperatures, leading to a Curie contribution to the susceptibility at small T. In 
this regime, the results are rather insensitive to U and Ef (U = lOf here). Fig. 1(b) is 
the n — 2.1 case. We chose U = t* here, and the results have a stronger dependence on 
Ef, showing a downturn at moderate T, before the Curie-law divergence sets in for T — > 0. 
Neither of these cases have a charge-transfer transition, rather they always display a Curie- 
like divergence as T — > 0. As n is decreased to 2 and beyond, it is no longer necessary 
for there to be any / electrons remaining at T = 0, and x^ can vanish in that limit. In 
Fig. 1(c) we plot for n = 2 and U — 2t*, showing an asymmetric peak that increases in 
size and sharpens as the turnover temperature Ty decreases, and which still has a weaker 
Curie upturn at the lowest temperatures, because rif ^ as T — > 0. Here, the results depend 
on Ef and U, with the peak lowering in magnitude and broadening as U increases, and the 
low-temperature upturn becoming more prominent. The case n = 1.9 is shown in Fig. 1(d) 
for U = lOt*. Once again we have a dependence on Ef and U, but the low-temperature 
upturn has disappeared. In some cases x^ has a double-hump structure. Finally the low- 
density regime is plotted in Fig. 1(e) (n = 0.5, U = lOf ). Here, as in the high-density 
limit, the results are insensitive to U. The reason why the results do not depend sensitively 
on U for both the high and low-density limits is that when n/ is small, then the Falicov- 
Kimball model looks like a noninteracting free band, because there are no localized electrons 
for the conduction electrons to scatter off of. Similarly, in the high density regime, as nf 
approaches 1, the model also becomes noninteracting, since there is an /-electron at every 
site, and the conduction band is simply shifted in energy by U , which does not affect any of 
the qualitative physics. 

The model also exhibits a metamagnetic transition, for large enough magnetic fields. This 
is shown in Figure 2 for one case: n = 2, Ef = t* , U = It* , g = 4.5, and eight temperatures 
from 0.05t* up to 6At*. The numerical analysis shows that at T = 0, the Zeeman energy 
at the transition is of the order of ksTy. The metamagnetic field, Hy, decreases at higher 
temperatures but the transition becomes smoother. The metamagnetism disappears as T is 
increased beyond Ty, because the /-electrons are already occupied in the zero-field limit. 

This anomalous behavior can be understood in terms of simple thermodynamic consid- 
erations. For Ef > ix and n < 2, the nonmagnetic empty /-orbital is energetically more 
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favorable than the magnetic state, so the ground state has no /-electrons. However, at 
high T, the large magnetic entropy of the /-electron spins favors the magnetic state. The 
susceptibility of these /-electrons grows rapidly as the temperature is reduced until, close 
to Ty, the entropy gain is insufficient to compensate the energy loss and the /-electrons 
disappear to form a nonmagnetic state. Thus, x drops rapidly, which is the origin of the 
sharp asymmetric peak in the x{T) data. 

IV. APPLICATION TO REAL MATERIALS 

The insulating state of N1I2 is best described by a hole picture. The localized holes 
lie in the Ni <i-bands, and there is one hole per site. The Iodine p-band is full, so there 
are no conduction holes. Hence the total hole concentration is 1 and in the regime where 
the pressure is tuned to lie just above the metallization pressure at T = 0, we expect the 
susceptibility to look qualitatively like the low-density regime of Fig. 1(e) (this is because the 
only difference between the hole picture and the electron picture is that the /-level energy 
changes sign). So, we expect an interesting anomalous magnetic response of MI2 when 
the pressure is sufficiently high that the system is just on the metallic side of the metal- 
insulator transition point. In this case, there are few localized holes at low temperature, 
but they will be thermally excited in large numbers as the temperature is raised. The 
susceptibility should be quite small at low temperature, and then rise dramatically as the 
metal-insulator transition point is passed (before turning over again at higher temperature). 
Also, we predict that MI2 should exhibit metamagnetism. Since new techniques allow one to 
inductively measure the magnetic response of samples in a diamond-anvil cell,0 without the 
need of attaching any wires within the cell, it is possible that one could use this anomalous 
magnetic response to accurately measure the phase diagram of N1I2. To our knowledge no 
measurements of the magnetic response of MI2 under pressure have yet been carried out. 
Only the Neel temperature in the insulating phase was measured with the Mossbauer effect.il 

For the Yb-based VF compounds, the mapping of this system to an effective Falicov- 
Kimball model is more complicated than in Nil 2 , and more controversial, so we begin by 
motivating how such a mapping can be made, and provide experimental evidence that 
supports this novel point of view. We notice first that there is no simple relation between the 
change in the lattice parameter (that is, the change in Yb-In hybridization) and the onset 
of the anomalies. In Ybi_ rr Y I InCu4 an increase of a(x) is accompanied by a decrease of Ty 
and Hy. In Ybln^^Ag^Cu^ a(x) does not change when In ions are replaced by smaller Ag 
ions (for x < 0.4) while Ty and Hy increase and the anomalies become less sharpHi. In 
stoichiometric YbInCu4, Ty and Hy are enhanced by thermal treatment, which introduces 
disorder between Yb and In sites but does not change a (or the Yb-In hybridization). We 
notice also that systems which are similar for T ^> Ty (H 3> Hy) and T <C Ty (H <C Hy) 
can show substantial variation in the onset and the shape of the anomalies at Ty and Hy. 
To account for such a behavior we assume that the effective f-state is in the proximity of 
the chemical potential, fi, so that correctly placed Yb ions are in the VF state. Hence, the 
important parameter is the renormalized position of the f-level, Ej — fi, which provides the 
characteristic energy scale (typically, about 100 K). 

The assumption that the majority of Yb ions are in the VF state for T 3> Ty is supported 
by the unusual deviation from Vegard's law found in YbIn 1 _ x Ag x Cu4 for x < 0.4 at room 



10 



temperature.! Here, the replacement of In ions by smaller Ag ions does not modify a(x), 
since in a VF compound the Yb ions reduce their average valence and preserve the volume of 
the unit cell. The assumption that Ef — // increases only slightly with x, explains the weak 
concentration dependence of the high-temperature susceptibility and the high-field satura- 
tion magnetization@i. For x < 0.4, Ef is pushed off which results in a HF compound 
with a stable f-shell.lll'! 

To explain the anomalies at Ty, we assume that the magnetic response of these Yb com- 
pounds has its origin in two physically distinct components. The first one, Xvf, is due to 
the hybridized (correctly placed) Yb ions, and is large (on an absolute scale), isotropic, and 
sample independent. The second component of the susceptibility, \ > arises from the unhy- 
bridized Yb ions which switch at Ty (or Hy) from the magnetic 3+ to the nonmagnetic 2+ 
configuration. Hence, x vanishes below Ty and is smaller than xvf at high temperatures, 
but dominates x(T) near T v ; it is strongly sample dependent. At low temperatures and 
high fields these unhybridized states lead to a metamagnetic transition with a small energy 
difference between the low-field and the high-field states. 

The lattice of YblnCui unit cells separates into two distinct sublattices: the sublattice 
of unhybridized Yb ions (localized f-electrons) and the sublattice of hybridized Yb ions. 
Conduction electrons from the In ions move (via nearest-neighbor hopping between unit 
cells of the full lattice) between any two unit cells of the unhybridized- Yb-ion sublattice. 
Hence, if we focus only on the sublattice degrees of freedom, we can model the d-band by 
a single effective band that couples any two sites of the sublattice with a random hopping 
matrix element (that depends on the growth and thermal treatment of the sample). These 
/ and d states have a common chemical potential (//) and interact by a Coulomb repulsion 
(U) when they occupy the same unit cell. In addition, both the d and / particles carry a 
spin label a and the d-level can accommodate 2 electrons (or holes) of the opposite spin. 
The occupancy of the /-level is restricted to rif < 1 because of the large Coulomb repulsion 
(Uff ~ oo) of the /-particles of opposite spins. To discuss the Yb-based VF compounds we 
use the hole-picture, in which Ef < \x and the total number of holes at the ill-placed sites 
is restricted to n h = n h d + n'j < 3. (In the electron-picture, one would have Ef > fj, and 
would restrict the total number of electrons to n e = n d + n'j < 3.) The magnetic field h 
couples to the / and d states but with different ^-factors (g = 4.5 for the /-holes). This 
picture is thus described by the spin one-half Falicov-Kimball model. In the limit N — > oo0, 
and N/Ni <C 1, our choice for maps this problem onto the infinite-dimensional (local) 
one, which allows the magnetic susceptibility to be evaluated exactly using the methods 
described above. 

The hole filling for these Yb-based VF compounds lies near n ~ 2. The exact solution 
reproduces well the overall behavior of the experimental data. Figs.(lc)-(ld) capture most 
of the features shownfHI! by x(T), such as an asymmetric peak at Ty that increases in 
magnitude and sharpens as Ty — > 0, while Fig. (2) explains the magnetizationi'0 M(T), 
with a metamagnetic transition that smoothes out and then disappears as T is increased. 
Various samples will have different numbers of Yb impurities and will require different values 
of t*, Ef, and n h . Our analysis shows that the lower the transition temperature, the more 
pronounced and steeper the anomaly. For large fields, unhybridized Yb ions switch to a 
magnetic configuration at much lower temperatures, than in the absence of the field. 

We emphasize that in the model proposed here, most of the Yb ions are in a hybridized 
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(VF) state at all temperatures and that the anomalies around Ty are due to an entropy 
driven transition of a small number of unhybridized Yb ions, so that the bulk of the lattice 
is unchanged at Ty. On the other hand, in the commonly used hybridization model it is 
assumed that all the Yb ions switch at Ty from a stable 3+ configuration to a hybridized VF 
configuration, which makes the high-T and the low-T phases fundamentally different. Thus, 
the appropriate model for describing the Yb-based VF compounds can be determined by 
measuring the pressure dependence of Aa/a across Ty. The Falicov-Kimball model predicts a 
weak pressure dependence because the ions involved in the transition are unhybridized, while 
the hybridization model predicts a strong pressure dependence because the hybridization 
rapidly increases with pressure. 

V. CONCLUSIONS 

In summary, we have exactly solved the spin-one-half Falicov-Kimball model on an 
infinite-dimensional hypercubic lattice in an external magnetic field. We also examined 
the magnetic response (as a function of temperature) in zero magnetic field. The system 
showed anomalous behavior due to the proximity of the /-electron states (which have local 
moments) to the chemical potential, which allows their occupancy to change dramatically 
as the temperature changes, because of the entropy gain due to the magnetic moments. 
Thus, the uniform magnetic susceptibility shows an asymmetric peak at a characteristic 
temperature Ty, which decays like a Curie- Weiss law for high temperatures (because of the 
/-moments) and typically decays exponentially rapidly for temperatures much lower than 
Ty (because the /-electrons are thermally excited across a gap). The magnetization (in 
an external magnetic field) shows metamagnetic behavior because of the rapid switching of 
electrons from the d- to the /-states when T < Ty, but shows a much smoother increase 
with magnetic field for temperatures larger than Ty. 

We applied this model to two candidate real materials. The first is Nil 2 , which is known 
to undergo a charge-transfer metal-insulator transition as a function of pressure, that is 
well described by the Falicov-Kimball model. We propose that the metal-insulator phase 
diagram, and this anomalous magnetic behavior can be easily measured using newly devel- 
oped techniques within a diamond-anvil cell. The second are the Yb-based VF compounds 
whose anomalous magnetic properties can be attributed to an entropy-driven transition of 
disordered Yb and explained by the Falicov-Kimball model with random hopping. The 
Falicov-Kimball-transition of unhybridized Yb ions and the VF behavior of hybridized Yb 
ions are both due to the proximity of the f-level to the chemical potential. The two inequiv- 
alent Yb sites appear in these compounds because of their characteristic crystal structure. 
Thus, the large changes in magnetic, transport, and elastic properties at Ty and Hy are 
explained, and reconciled with small changes in the /-electron valence and Aa. This pro- 
posed alternative explanation for the magnetic response of the Yb-based VF compounds is 
controversial, and differs from the more conventional description (which does not have any 
quantitative theoretical description). We propose high-pressure experiments to differentiate 
between the two proposed theories. 
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FIGURES 



FIG. 1. Magnetic susceptibility of the /-electrons in the Falicov-Kimball model. The different 
curves correspond to Ef/t* = 0.0,0.5, . . . ,4.5 (in general Ef increases from top to bottom in these 
figures). The different figures are: (a) n = 2.5 and U = Wt*; (b) n = 2.1 and U = t*; (c) n = 2.0 
and U = 2t*; (d) n = 1.9 and U = Wt*; and (e) n = 0.5 and U = Wt*. 

FIG. 2. Total magnetization of the Falicov-Kimball model as a function of magnetic field. 
The different curves correspond to different temperatures: T/t* = 0.05,0.1,0.2,... ,3.2,6.4 (the 
temperature increases from top to bottom in the large- ft, range). The parameters are n = 2, 
E f = t*, and U = 2t*. 
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